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(57) Abstract: The present invention is related to a method For performing digital holographic imaging (DIN), said DHI being 
characterized by the fact that images of a sample (0 are obtained by applying numerical means to reconstruct holograms of the 
sample. In the method of the invention, the sample is in a medium (3) with controlled properties (5) that influence the behaviour of 
the sample and/or that influence the process of hologram formation. The information content of one hologram, or of a plurality of 
holograms' recorded with the sample in one medium, or in a plurality of medium, is used to reconstruct one image ol the sample, or 
a plurality of images of the sample, that describe quantitatively one property of the sample, or a plurality of properties of the sample. 
The present invention is also related to an apparatus with which to perform said method. 
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Apparatus and method for digital holographic imaging 



Field of the invention 

5 The invention relates to digital holographic imaging (DHI) apparatuses and methods which 
provide an holographic representation of a sample. 

Sate of the art 

Holography is a three-dimensional (3D) imaging technique that makes use of the 

10 interference between a reference wave and a wave emanating from the sample called object 
wave. The purpose of this interference is to record the phase of the object wave, which is 
related to the 3D character of the sample. With digital holographic imaging (DHI), real-time 
observations can be achieved by using a charged coupled device (CCD) camera as recording 
device and by performing a numerical reconstruction of the hologram. This idea has been 

15 proposed for the first time over 30 years ago by J. W. Goodmarm, R. W. Lawrence, in 
"Digital image formation from electronically detected holograms," Appl. Phys. Lett, Vol. 11, 
1967. As a result of technological progresses achieved in the fields of digital image 
acquisition and processing, this numerical or digital approach of holography has considerably 
extended the fields of its potential applications and different types of DHI-inspired imaging 

20 systems have been developed during the last years. 

DHI techniques can be classified in two main categories: in-line techniques 
characterized by the fact that the reference and object waves have similar propagation 
directions, and off-axis techniques for which the two interfering waves propagates along 
different direction. The procedure for hologram formation in in-line digital holography is 

25 similar to the procedure used for phase measurements with so-called phase-shifting 

1 
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interferometric techniques. Hologram formation with in-line techniques requires the 
acquisition of several images, at least three, that must be recorded during a modulation of the 
reference phase. Off-axis techniques, are more simple from the experimental point of view 
since they require a single hologram acquisition without modulation of the phase of the 

5 reference wave. In-line techniques however present the advantage that the reconstructed 
images are free of twin images and zero order of diffraction. Among off-axis techniques, we 
can distinguish methods based on Fourier-transform holography, and methods based on a so- 
called Fresnel holography. With Fourier-transform methods the reference wave must be a 
spherical wave of precisely controlled curvature and image reconstruction is basically 

10 performed by Fourier transformation of the hologram. With Fresnel-holography based 
techniques, the reconstruction procedure is more sophisticated but more flexibility is offered 
to build experimental installations. 

Among recent publications presenting developments or applications of DHI-inspired 
techniques, we can mention the following works. A study of some general performances of an 

15 in-line technique is presented in "Image formation in phase-shifting digital holography and 
application to microscopy", I. Yamaguchi et al., Applied Optics, Vol. 40, No. 34, 2001, pp. 
6177-6186. In "Fourier-transform holographic microscope," Applied Optics, Vol 31, 1992, 
pp. 4973-4978, W.S. Haddad et al describe the general principle of Fourier-transform DHL 
Examples of applications of the Fresnel-based approach can be found in "Direct recording of 

20 holograms by a CCD target and numerical reconstruction," U. Schnars and W. Juptner, 
Applied Optics, Vol. 33, 1994, pp. 179-181, and in "Performances of endoscopic holography 
with a multicore optical fiber," O. Coquoz et al., Applied Optics, Vol. 34, 1995, pp. 7186- 
7193. 

A key element of a DHI method is the numerical method used for hologram 
25 reconstruction. An original reconstruction procedure, which allows for reconstructing 
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simultaneously the amplitude and tne pnase or me ooject wave, on inc oasis 01 a single on- 
axis hologram acquisition, has been developed by Cuche et al and is presented in US patent 
6,262,218, and in WO 00/20929. Different applications and implementations of this technique 
are presented in "Digital holography for quantitative phase-contrast imaging", Optics Letters, 

5 Vol. 24, 1999, pp. 291-293, in "Simultaneous amplitude-contrast and quantitative phase- 
contrast microscopy by numerical reconstruction of Fresnel off-axis holograms", Applied 
Optics, Vol. 38, 1999, pp. 6994-7001, in "Spatial Filtering for Zero-Order and Twin-Image 
Elimination in Digital Off-Axis Holography", Applied Optics, Vol. 38 No. 34, 1999, in 
"Aperture apodization using cubic spline interpolation: Application in digital holographic 

10 microscopy", Optics Communications, Vol. 182, 2000, pp. 59-69, and in "Polarization 
Imaging by Use of Digital Holography", T. Colomb et al., Applied Optics, Vol. 38, No 34, 
1999. 

DHl method presents interesting possibilities of applications in cell biology. Indeed a 
living cell behaves optically as a phase object, i.e. a transparent sample whose constituents 

15 can be optically probed on the basis of the phase shift they induce on the light crossing them. 
The phase-shifting behavior of transparent sample is well known, and for a long time as it 
constitutes the mechanism of image formation in phase-contrast (PhC) and Nomarski (DIC) 
microscopy. Even though these two techniques are widely used in biological microscopy, and 
well suited as contrasting methods, they cannot be used for precise quantitative phase 

20 measurements. The DHI method instead, is reminiscent of classical interferometry, which is 
the most commonly used technique for phase measurements. However, whereas 
interferometric techniques are widely used in metrology, only few biological applications 
have been reported, by R. Barer and S. Joseph, in "Refractometry of living cells", Quarterly 
Journal of Microscopical Science, Vol. 95, 1954, pp. 399-423, by .R. Barer in "Refractometry 

25 and interferometry of living cells," Journal of the Optical Society of America, Vol. 47, 1957, 

3 
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pp. 545-556, by A. J. Coble et al. in "Microscope interferometry of necturus gallblader 
epithelium", Josiah Macy Jr. Fundation, New York, 1982, p. 270-303, by K. C. Svoboda et al 
in "Direct observation of kinesin stepping by optical trapping interferometry ", Nature, Vol. 
365, 1993, by J. Farinas and A. S. Verkman, in "Cell volume plasma membrane osmotic 
5 water permeability in epithelial cell layers measured by interferometry," Biophysical Journal, 
Vol. 71, 1996, pp. 3511-3522, by G.A. Dunn and D. Zicha in "Dynamics of fibroblast 
spreading," Journal of Cell Science, Vol. 108, 1995, pp. 1239-1249. 

For biological applications, as well as for material science or metrology applications, 
DHI methods offer a novel alternative to classical interferometry with similar performances 

10 but simplified experimental procedures. The main advantage originates from the fact that 
complex and costly experimental optical devices can be handled by digital processing 
methods. For example, as described by E. Cuche et al. in "Simultaneous amplitude-contrast 
and quantitative phase-contrast microscopy by numerical reconstruction of Fresnel off-axis 
holograms", Applied Optics, Vol. 38, 1999, pp. 6994-7001, the wave front deformations 

15 appearing when a microscope objective is introduced along the path of the object wave can be 
compensated using a digital procedure. This particular feature opens attractive possibilities in 
the fields of microscopy. In addition DHI techniques performs faster than interferomeric 
techniques, and provides more information about the sample, in particular, the amplitude and 
the phase of the object wave can be obtained simultaneously on the basis of a single hologram 

20 acquisition. 

DHI methods have been applied to static imaging of biological cells, without phase 
reconstruction by K. Boyer et al. in "Biomedical three-dimensional holographic microimaging 
at visible, ultraviolet and X-ray wavelength", Nature Medecine, Vol. 2, 1996, pp. 939-941, 
and by F. Dubois et al. in "Improved three-dimensional imaging with a digital holography 
25 microscope with a source of partial spatial coherence," Applied Optics, Vol. 38, 1999, pp. 
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7US!>-70y4. DHI ot cells using a pnase measurement moaauty requinng several image 
acquisitions has been reported by G. Indebetouw and P. Klysubun in "Saptiotemporal digital 
microholography," Journal of the Optical Society of America A, Vol. 18, 2001, pp. 319-325. 

With DHI, image acquisition can be performed at video-rate, and even faster using 
appropriate image acquisition systems, for experimental periods of up to several hours. Due to 
its interferometric nature, DHI has a high axial resolution (nanometer scale), which allows for 
observing subtle and minute modifications of sample shape, opening a wide field of 
applications in both life and material sciences. 

A well-known limitation of interferometric techniques in general, and of DHI methods 
in particular, comes from the fact that the measured quantity, i.e. the phase of an 
electromagnetic radiation, can be simultaneously influenced by several factors. When an 
electromagnetic radiation interact with a sample, the result of this interaction, and in particular 
the changes occurring on the phase of the incident radiation, depends on three main classes of 
parameters: 

a) Geometrical properties of the sample . The surface topography for a reflective 
sample imaged in reflection and the thickness of the sample for a transparent 
sample imaged in transmission. 

b) Dielectric properties of the sample. The absorption that influences the intensity 
of the radiation, and the index of refraction that modifies the phase of the 
radiation. Absorption and index of refraction are often combined in the concept 
of complex index of refraction. 

c) Dielectric properties of the medium incorporating the sample. It is only in the 
vacuum, or by approximation in the air, that an electromagnetic interaction 
depends only on the properties of the sample. In the general case, the dielectric 
properties of a sample must be defined relatively to the dielectric properties of 

5 
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the medium incorporating the sample. For example, its is not the absolute 
index of refraction of the sample (n s ) that must be considered, but the 
difference (n s -n m ) f where n m is the index of refraction of the medium 
incorporating the sample. 
5 It is therefore often difficult to interpret data resulting from the measurement of the phase of 
an electromagnetic radiation. For example, if one wants to determine precisely and 
unambiguously the shape of a sample, the knowledge of the dielectric properties (of both the 
sample and the medium) is required. And if one wants to determine the dielectric properties of 
a sample, the shape of the sample and the dielectric properties of the medium must be known. 
10 A rigorous and unambiguous interpretation of the results can be obtained for example by 
taking dielectric property values from the literature, or by measuring the geometry or the 
dielectric properties of the sample with another technique. 

Objectives of the invention 

15 The present invention aims at providing an apparatus and a method for performing 

digital holographic imaging (DHI) which consider the medium incorporating the sample as an 
active parameter of the image formation process. It is also an objective that such an apparatus 
and method improve the performances of DHI by providing more accurate measurements. It 
is another objective to study dynamic effects resulting from interactions between the sample 

20 and the incorporating medium. Finally, another objective is to provide more information 
about the sample, in particular concerning its geometrical and dielectric properties. 



25 
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Summary 01 the invention 

The present invention concerns an apparatus as defined in claim 1, a medium as defined 
in claim 6 and a method for performing digital holographic imaging (DHI) of a sample as 
defined in claim 7. 

5 Some embodiments of the invention which relate to an apparatus are defined in 

subclaims 2 to 5. Other embodiments which relate to a method are defined in subclaim 8 to 
13. 

The invention is particularly useful with a method or a device as disclosed in PCT patent 
application WO 00/20929 and in particular with a device comprising an hologram creation 

10 device which includes a source of radiation, means to illuminate the sample unit with said 
radiation, means to create an object wave by collecting the radiation after interaction with the 
sample unit, means to create a reference wave from said radiation, means to produce a 
hologram by interference between said reference wave and said object wave, means to 
produce a digital hologram by acquiring and digitizing said hologram, and means to transmit 

15 said digital hologram to said hologram reconstruction unit, and wherein said hologram 
reconstruction unit provides an amplitude-contrast image and/or a quantitative phase contrast 
image by numerical calculations applied to said digital hologram, and wherein said processing 
unit determine the dielectric properties of the sample and/or the shape of the sample by 
processing a plurality of said amplitude contrast images and/or a plurality of said quantitative 

20 phase-contrast images, and wherein said processing unit can also be used to describe the 
evolution of said sample by processing a plurality of said amplitude contrast images and/or a 
plurality of said quantitative phase-contrast images. 

25 
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Short description of the drawings 

Fig. 1 is a schematic diagram showing the building block of an apparatus according to the 
present invention. 

Fig. 2 refractive index (n) and absorption coefficient (k) behaviors near an absorption band. 

5 Fig. 3 presents an example of application of the method of the present invention. Fig 3a 
schematically represent a biological cell in a perfusion chamber. Fig 3b shows an example of 
a hologram of a biological cell. Fig. 3c shows an example of a quantitative phase-contrast 
image of a biological cell. Fig. 3d shows an example of a quantitative phase-contrast image of 
a biological cell presented in 3D perspective. 

10 Fig. 4 presents an example of application of the method of the present invention. Fig 4a shows 
an example of the evolution during time of the phase measured in a specific area of a cell. Fig. 
4b shows an example of a quantitative phase-contrast image of a cell. Fig 4c shows an 
example of a concentration-response curve established thanks to the present invention. 

15 Detailed description of the invention 

An important aspect of the invention is represented by the medium which includes the 
sample. As schematically shown in Fig. 1, this means that the sample 1 is located inside a 
given volume, defined for example by a container 2, filled with a medium 3. The sample unit 
4 is composed of the sample 1, the container 2 and the medium 3. The medium 3 can be a gas, 
20 a liquid, a gel, a solid, a powder or a mixture of solid particles in suspension in a liquid. If the 
medium 3 is a gas, the container 2 is closed. If medium 3 is a liquid, the container 2 can also 
be closed, but the top face of the container 2 can be opened. If medium 3 is solid or a gel, 
container 2 is optional. For example, the sample 1 can be a living sample embedded in a 
physiological solution. 
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the index of refraction of medium 3. For example, medium 3 can be a liquid element 
comprising different concentrations of molecules such as metrizamide (C18H22I3N3O8), and/or 
different concentrations of manitol, and/or iohexol, and/or iodixanol, and/or ficoll, and/or 
5 percoll. 

In particular, medium 3 may be composed of, or may comprise, elements that have 
dispersion properties, meaning that the dielectric properties of medium 3 are different for 
different wavelength of the electromagnetic spectrum. For example, medium 3 may comprise 
concentrations of an organic molecule, and/or of an inorganic molecule, and/or of a dye, 

10 and/or of a stain, and/or of a chromophore, and/or of pigments, and/or of a colloid, and/or of a 
neutral stain, and/or of a diazotic compound, and/or of a triphenylmethan compound, and/or 
of a xanthen compound, and/or of a aminoxanthen compound, and/or of a hydroxyxanthen 
compound, and/or of pyronyn, and/or of rhodamin, and/or of fluorescein, and/or of eosin, 
and/or of carbocyanin, and/or of oxadicarbocyanin, and/or of methylen blue, and/or of 

15 phthalocyanin, and/or of a ph-indicator element, and/or of sulforhodamin, and/or of trypan 
blue, and/or of thiazin, and/or of fast green FCF, and/or of evan's blue. 

In particular, medium 3 may be a mixture of particles, micro- and/or nano-particles in 
particular, in suspension in a liquid. For example, these particles may be metallic particles, 
and/or dielectric particles, and/or coated particles, and/or gold particles, and/or silver 

20 particles, and/or copper particles, and/or silicon particles, and/or crystalline particles, and/or 
semi-conductors particles, and/or CdS particles, and/or CdSe particles, and/or latex beds, 
and/or charged particles, and/or colloidal gold particles, and/or silver particles. 

A medium controller device 5 can be optionally used, more particularly if medium 3 is 
a gas or a liquid. The role of the medium controller device 5 is to permit exchanges of 

25 medium 3. When medium 3 is composed of a plurality of constituents, the medium controller 
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device 5 can be used to change the concentration of one of these constituents, or of a plurality 
of these constituents. The medium controller device 5 can also comprise means to control the 
temperature, and/or the pressure of medium 3. It may also include means to control gas partial 
pressures such as oxygen or carbon dioxide. The medium controller device 5 can also 
comprise means, to induce electrical or magnetic stimulations of the sample 1 or of the entire 
sample unit 4. In particular, electrodes or electrode arrays may be included in the medium 
control device 5. The medium controller device 5 can also comprise means to induce a 
mechanical stimulation of the sample, for example a tip and a cantilever as performed in 
atomic force microscopy (AFM). 

In particular, medium 3 may be a mixture of particles, micro- and/or nano-particles in 
particular, in suspension in a liquid, and the medium controller device 5 may be used to 
control the type of these particles, and/or the concentration of these particles, and/or the sizes 
of these particles. 

When introduced in a set-up for digital holographic imaging (DHI), the sample unit 4 
is illuminated by an incident radiation called illumination wave 6. After interaction with the 
sample unit, the incident radiation becomes the object wave 7. The incident radiation is 
provided by a radiation source 8, which emits preferentially radiations of the electromagnetic 
type. It is clear that to enable interaction between the sample 1 and the illumination wave 6, at 
least one face of container 2 comprises a zone that is transparent for the incident radiation. 

The medium controller device 5 may comprise means, such as piezo-electric or 
acoustic transducers, to induce mechanical vibrations of the sample unit 4 or a vibration of 
medium 3. Such vibrations are expected to have a positive impact on the signal to noise ratio 
(SNR) of the method by inducing a random variation of the propagation directions of parasitic 
radiations emitted by elements of the container 2 that interact with the illumination wave 6. 
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example in patent WO 00/20929. We just give here a summary of the most important steps of 
the DHI process. The object wave 7 is collected using appropriate means 9 (e.g. a lens, an 
assembly of lenses or a microscope objective), and appropriate means (e.g. an assembly beam 

5 splitters, mirrors, lenses and apertures) are used to produce a reference wave 10 from the 
radiation source 8. Appropriate means 1 1 (e.g. an assembly beam splitters, mirrors, lenses and 
apertures) are used to form a hologram 12 by making the reference wave 10, and the object 
wave 7., to interfere with each other. The hologram 12 is acquired by an image acquisition 
system 13, which transmit a digital hologram 14 to a hologram reconstruction unit 15. The 

10 reconstruction unit 15 applies numerical calculations to reconstruct an amplitude contrast 
image 16 and/or a quantitative phase contrast image 17. An additional processing unit 18 can 
be used to collect and to process a plurality of said amplitude contrast images and/or a 
plurality of said quantitative phase contrast images. The main function of the processing unit 
18 is to provide results that are specific to the method of the present invention. In particular, 

15 the processing unit 18 can be used to describe the evolution of the sample 1 during time, 
and/or to provide measurements and a representation its dielectric properties, and/or to 
provide measurements and a representation of its thickness, and/or to provide measurements 
and a representation of its topography. 

The processing unit 18 can also be used to improve the accuracy of the measurements 

20 by applying statistical procedure to a plurality of said amplitude contrast images 16 and/or to 
a plurality of said quantitative phase contrast images 17. 

When a plurality of amplitude contrast images 16 and/or a plurality of quantitative 
phase contrast images 17 are provided during a given period of time, the processing unit 18 
can be used to process temporal signals extracted from the plurality of amplitude contrast 

25 images 16 and/or the plurality of quantitative phase contrast images. For example, the 

11 
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processing unit can be used to Fourier-transform said temporal signals, and/or can be used to 
apply band-pass, and/or low-pass, and/or high pass filtering techniques to said temporal 
signal. 

Alternatively a plurality of amplitude contrast images 16 and/or the plurality of 

5 quantitative phase contrast images 17 can be collected after successive modifications of the 
irradiation wavelength X. If the holograms are recorded in the transmission geometry with a 
semi transparent sample, and provided that the dispersion law of the refractive index is 
known, a tomographic image of the sample can be obtained 

In another alternative a plurality of amplitude contrast images 16 and/or the plurality 

10 of quantitative phase contrast images 17 can be collected at variable incidence angles of the 
irradiating beam. This achievement can be obtained by rotating the sample. In another 
embodiment, the sample can be kept fixed and the irradiating beam can be moved in order to 
vary the incidence angle of the irradiating beam over the sample. The processing unit 18 can 
be used to operate a Radon transform of the collected data. Furthermore, such a mean can be 

15 operated to compute the local dielectric properties of the sample by applying an algorithm of 
back-projection. Fourier slice theorem in particular can be applied to the calculation of the 
dielectric potential of the sample. 

In another alternative a plurality of amplitude contrast images 16 and/or the plurality 
of quantitative phase contrast images 17 can be collected after successive modifications of the 

20 irradiation wavelength X and after successive modifications of the incidence angles of the 
irradiating beam. The processing unit 18 can be used to compute the local dielectric properties 
of the sample by applying an algorithm. Fourier slice theorem as commonly applied in 
diffraction tomography. 

The processing unit 18 can also be used to generate representations of the obtained 

25 results, for example by producing three-dimensional perspective representations of the 
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images, and/or by generating movies or animations from a plurality of obtained images. 

When a plurality of images are provided, the processing control unit can also be used 
to compare them which each other, for example in order to reveal the spatial distribution of 
5 signal changes. 

The most important physical parameters describing a typical application of the present 
invention are : 

- § is the phase of the object wave 7, expressed in degree or in gradient. <|> is provided 
by DHI as an image called quantitative phase-contrast image 17. 

10 - OPL=§)J2n , where X is the wavelength of the radiation. OPL is called optical path 

length. When <j> is measured in degree, we have OPL= 

- n m is the index of refraction, also called refractive index, of the medium 
incorporating the sample, or more precisely the index of refraction of medium 3. 

- n s the index of refraction of the sample 1 . 
15 - d the thickness of the sample 1 . 

- /; the height of the sample 1 , which describes the topography of the sample. 

All parameters described in the list above are provided in the form two-dimensional (2D) 
functions describing the spatial distribution of these parameters in a specific plane, or in other 
words in the form of images. It is clear when these images are provided in a digital form, the 
20 corresponding parameter can be measured from the values of the pixels composing the digital 
images. 

In a transmission imaging geometry, which is characterized by the fact that the object 
wave 7 is collected after the illumination wave 6 has crossed the sample unit 7, the quantity 
measured in one pixel of said quantitative phase-contrast image maps, the 2D distribution of 
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the phase shifts induced by the samplel on the radiation crossing the sample 1. The 

corresponding OPL can be expressed as follows: 

OPL = (n 5 -n m )d (1), 
Expression (1) clearly indicates that the measured signal depends on the structural properties 
5 of the sample (its thickness d), and on the difference between the dielectric properties of the 

sample, and the dielectric properties of the medium incorporating the sample (n s - n m the 
difference of index of refraction). 

In a reflection imaging geometry, characterized by the fact that the object wave 7 is 
collected after the illumination wave 6 has been reflected or backscattered by the sample, and 

10 if we consider only sample having surfaces with constant and homogeneous dielectric 
properties, the quantity measured in one pixel of said quantitative phase-contrast image 
depends on the surface topography of the sample and on the index of refraction of the medium 
incorporating the sample n m . The corresponding OPL can be expressed as follows: 

OPL = 2n m d (2), 

15 Expressions (I) and (2) gives only examples of expressions describing the behavior of 

OPL as a function of the specimen shape and dielectric properties. The present invention can 
also be applied determine the specimen shape and the dielectric properties of the specimen on 
the bases of more complex or more general models describing the interaction of an 
electromagnetic radiation with a specimen. In particular, such models can be obtained within 

20 the framework of the Maxwell equations, and/or within the framework of a vectorial or 
tensorial theory of electromagnetic interaction, and/or within the framework of the scattering 
potential equation, and/or within the framework of theoretical model describing the 
diffraction of an electromagnetic radiation by a specimen. It is also clear that such models 
may involve other dielectric properties than the refractive index, in particular the absorption, 

25 and/or the complex index of refraction. In particular such models may give more accurate 
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dielectric properties of medium are significantly different. 

It is an object of the present invention that the proposed method and apparatus can be 
used to study the effect of a medium on the sample. In this case, the medium controller device 

5 5 can be used to control the quantity of a given medium in the medium incorporating the 
sample. If this medium interact with the sample and induces OPL variations, theses OPL 
variations can be measured and monitored during time with DHL 

It is an object of the present invention that the medium incorporating the sample can 
be used to modulate the OPL signal measured by DHL More particularly, and as clearly stated 

10 by expressions (1) and (2), changing n m has a direct influence on the measured OPL signal. In 
particular, for small or thin sample producing weak OPL signals, an OPL signal amplification 
can be achieved by increasing //,„ in reflection, or by increasing n s - n m in transmission. 

Changing n m the index of refraction of the medium incorporating the sample can be 
achieved by changing medium 3 using the medium controller device 5. In particular, the 

15 medium controller device 5 can be used to control the composition of medium 3, said 
composition comprising elements that changes the index of refraction of the medium. For 
example, medium 3 can be a liquid element comprising different concentrations of molecules 
or ions. For example medium 3 can comprise elements such as metrizamide (C18H22I3N3O8), 
and/or manitol, and/or iohexol, and/or iodixanol, and/or ficoll, and/or percoll. For example, 

20 medium 3 can be a mixture of particles in suspension in a liquid, and the index of refraction 
of the medium 3 can be changed by changing the distribution of the sizes of these particles, 
and/or by changing the statistical distribution of the size of these particles, and/or by changing 
the concentration of these particles, and/or by changing the material of these particles. 

It is an object of the present invention that the medium 3 enveloping the sample I may 

25 be composed of, or may comprise, an element having dispersion properties, meaning that the 

15 
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index of refraction of the medium incorporating the sample n m is different for different 
wavelength of the radiation. Therefore, n m can be changed by changing the wavelength X of 
the radiation illuminating the sample unit 4. 

For example the medium having dispersion properties can be characterized by an 

5 absorption spectrum showing absorption lines or absorption bands. According to the Kramers- 
Kronig relationship the dispersion behavior of such a medium, in a wavelength interval near 
an absorption line or near an absorption band, is characterized by the fact that important 
refractive index changes can occur even in small wavelength intervals. For example, as 
described in Fig. 2, two wavelengths X\ and *k 2y can be chosen so that they are located along 

10 the decreasing side of an absorption band, where dispersion effects are more pronounced. 
This enables fast and important modulations of n m . 

Particularly in the transmission imaging geometry, when a plurality of quantitative 
phase-contrast images have been obtained for different values of n m , the thickness of the 
sample (rf), and the index of refraction of the sample n Si can be obtained separately or 

15 simultaneously, assuming that the different values of n m are known or have been measured 
previously. For example, lets define OPL1 and OPL2, two different OPL measurements 
performed for two different values of n m , respectively n ml and n m2 . From expression (1), we 
have: 

20 OPL\-(n s -'n mi )d (3), 

OPL2 = (n s - n m2 )d (4), 

and n s and d can be directly determined as follows: 

25 d = (OPL1 - OPL2;/(/i m2 - n m i) (5), 

n s = (OPL2 n mi - OPL1 n „ l2 )/(OPL2 - OPL1) (6). 
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Jiven n iwo image acquismons are suiiiciern, mure auquisiuuns eau uc pcuuuucu ioi 
different values of n m , or for the same value of n mj for example in order to increase the 
accuracy of the measurements. 

Different values of n m (for example the two values n ml and n m2 ) can be obtained by 
5 changing the medium. Different values of n m (for example the two values n m! and n m2 ) can 
also be obtained by using a medium having dispersion properties that allows to change its 
index of refraction by changing the wavelength of the radiation illuminating the specimen 
unit. 

The invention can be used to measure «,„. In the transmission imaging geometry for 
10 example, if d and n Sf are known, n m can be determined as follows: n m = n 5 -OPL/d. More 
generally, the present invention can be used as a technique for refractometry when the 
geometry of the specimen is known. 

The invention can also be used to obtain a tomography of the sample. Provided that 
n m , can be chosen reasonably close to n s in order to avoid significant refraction effects, and 
15 provided that n m , and n s Aon not vary to much over a distance comparable to the wavelength A., 
in order to avoid significant diffraction effects, we have in a first approximation for a 2D 
geometry:: 

OPL®,r)= JJ dxdy n (x y y)b(x cos8 + y sin9 -r) (7) 

Specimen 
unit 

where n(x,y) is the local refractive index which is equal to n m outside the sample and 
20 approximately n 5 with local variations over the volume of the sample. 5 is the Dirac delta 
distribution. OPLfi.RJddXz can be obtained by collecting a plurality of amplitude contrast 
images 16 and/or a plurality of quantitative phase contrast images 17 can be collected at 
variable incidence angles 9 of the irradiating beam. From the collected data the refractive 
index n(x.y) of the sample can be computed by applying an algorithm of back-projection. 

17 
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If medium 3 may have no influences on the sample, the invention can be used to 
measure a property or a behavior of the medium 3. For example, if medium 3 is a liquid, the 
present invention can be used to study, or to observe, dynamics fluid phenomenon introducing 
local variations of OPL, flow analysis in particular, for example with flows tracers. For 

5 example, and more particularly when medium 3 is a solid or a gel, the present invention can 
also be used to analyze internal stresses inside the material. 

As described by T. Colomb et al. in "Polarization Imaging by Use of Digital 
Holography", Applied Optics, Vol. 38, No 34, 1999, DHI can be used to image and to 
measure the polarization sate of the object wave. It is an object of the present invention that, 

10 in addition to amplitude-contrast images 16, and quantitative phase-contrast images 17, the 
polarization sate of the object wave, may be provided to the processing unit 18. In this case, 
the present invention can also be used to observe and to analyses processes that influence the 
polarization of the radiation, and in particular dynamic processes inducing transient 
modifications of the polarization. In this case for example, the' medium 3 may have 

15 birefringent properties. 

If a plurality of holograms are acquired during a period of time, the corresponding 
plurality of quantitative phase contrast images can be processed, for example by the 
processing unit 18, in order to compensate temporal variation of the OPL or phase signal, said 
temporal variations having other origins than the behavior of the sample, for example such 

20 temporal variation may be caused by mechanical vibrations of elements composing the 
apparatus used to generate holograms, and/or may be caused by air turbulences and/or may be 
caused by thermal effects. To compensate such temporal variations, the present invention uses 
' a computer-based method that measures the phase signal in regions of the quantitative phase- 
contrast images where the sample has no influence on the OPL signal, or in other words in 

25 regions where d, or h, or n m have no influence on the measured signal. In such regions, OPL 
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measurements provide a reterence signal mai monuor me miiuence 01 unucsnauic sources 01 
signal variations. The compensation of the influence of these undesirable sources of signal 
variations can be performed by subtracting the reference signal from the signal measured in 
other area of the quantitative phase-contrast images, in particular in area where the sample 
influences the measured signal. 

Within the general context of imaging and microscopy techniques, DHI methods in 
general, and the present invention in particular, appear as a powerful mean for real-time 
imaging of processes inducing dynamic changes of a sample shape and/or dynamic changes 
of physical properties of a sample. These dynamic changes may appear in relation to a 
spontaneous process and/or in relation to an induced process. For example, DHI methods can 
be applied to monitor the response of a sample to external perturbations of different forms, for 
example chemical perturbations, and/or biological perturbations, and/or electrical 
perturbations, and/or mechanical perturbations. 

DHI methods in general, and the present invention in particular, present attractive 
possibilities for multifunctional imaging approaches combining DHI with other types of 
existing imaging techniques, and/or with other types of measurement techniques. For 
example, the present invention can be combined with fluorescence microscopy, or with laser 
scanning confocal microscopy, or with two-photon microscopy, or with second or third 
harmonic generation microscopy, or with atomic force microscopy. The present invention can 
also be applied in parallel with experimental techniques used in electrophysiology, and/or in 
fluid mechanics. 

DHI techniques in general, and the present invention in particular, are particularly well 
adapted for non-invasive analysis, because no contrast agents are required, and because low 
radiation intensities are used. In particular, it has been established that the measured 
irradiance at the sample plane may be several orders of magnitude lower than the irradiance 

19 
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used for example by other microscopy or imaging techniques. It is therefore an interesting 
potential application of the present invention that the method can be applied to image delicate 
sample, and biological sample in particular. For example, the sample can be a biological cell, 
a living biological cell, a culture of biological cell, a mono-layer of cultured biological cells, a 

5 tissue composed of biological cell, a biochip, a preparations of proteins. 

A particular aspect of the method according to the invention is characterized by the fact 
that the medium incorporating the sample can be changed during time and/or that properties 
of the medium incorporating the sample can be changed during time. Such changes of the 
medium are expected to induce modifications of the sample shape and/or modifications of the 

10 dielectric properties of the sample, and said modifications of the sample are monitored by 
DHL 

As a result of the high acquisition frequency of the DHI method, the method of the 
present invention is particularly well adapted to study the dynamics of fast sample-medium 
reactions. As DHI is also a non-invasive imaging technique working without contrast agents 
15 and with low radiation intensities, the method of the present invention is particularly well 
adapted to study the dynamic of reactions with minimized external perturbations, and during 
long time periods of hours or even days. This is of particular interest for applications in life 
sciences. 

Fig. 3 presents an example of application of the present invention to image living cells 
20 in culture enclosed in a perfusion chamber 30 filled with a physiological solution 37. As 
shown in Fig. 3a, the cell culture 31 is illuminated by an incident radiation 32 and a 
microscope objective 33 collects the transmitted radiation to form an object wave 34. A CCD 
camera 35 records the hologram. An example of hologram is presented in Fig. 3b. The 
hologram is produced by the interference with a reference wave 36. A numerical procedure is 
25 then applied to reconstruct the hologram. An example of quantitative phase-contrast image is 
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presented in rig. sc. ine ootainea pnasc-coiurasi image can oe preseiucu in pcispecuve, 
as presented for example in Fig. 4d. 

From the quantitative phase-contrast image, quantitative data about the cell morphology 
can be obtained. The refractive index, at the wavelength of 633nm, of the perfused solution n m 
5 = 1.3336 ± 0.0002, was measured with an Abbe refractometer. By assuming, in first 
approximation, a constant and homogeneous cellular refractive index n s = 1.365, one can 
estimate that an phase of 10 degree corresponds to a cellular thickness d = 0.6 ^m. Using this 
relation between phase and thickness, phase measurements extracted from specific pixels of 
the quantitative phase-contrast image presented Fig. 3c, yields a thickness of 1-3 jim for the 
10 neuronal processes and of approximately 8 }im for the cell body. Even with the assumption of 
a homogeneous distribution of intracellular refractive index (which can be discussed given a 
certain degree of heterogeneity of its constituents), these estimations of cell morphology give 
realistic values of typical cell dimensions as described in the literature for the mouse cortical 
neurons presented here. Using the same assumption for refractive index of the cell, and based 
15 on a phase measurement accuracy of 0.7 degree, that is dictated by experimental opto- 
mechanical noise, we obtain an axial (or vertical) resolution of 40 nanometers. 

As an example of application of the present invention, the present invention has been 
applied to study the effects of an excitatory neurotransmitter (glutamate) on the morphology 
of cultured mouse cortical neurons. The influence of glutamate on neuron morphology has 
20 been already described as resulting from excitotoxicity, for high concentrations (20-500 \\M) 
and long exposures of several minutes. Thanks to its unique features, the method of the 
present invention has made possible the observation of reversible morphological changes 
occuring for short applications (a few seconds), at physiological concentrations (EC 5 o = 7.5 
^M). It has also been established that this effect occurs in response to the activation of N- 
25 methyl-D-aspartate (NMDA) receptors. Results of the real-time observation of glutamate- 
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mediated morphological changes in neurons are summarized in Figs. 4 and 5. Figure 4a 
presents the temporal evolution of the phase signal measured, in real-time by DHI, in the 
domain delimited by a black rectangle in Fig. 4b, which presents a typical example of a 
quantitative phase contrast image of a living mouse cortical neuron immersed in a 
5 physiological solution. Cells were observed in a control solution (HEPES-buffered 

physiological solution) and pharmacological agents were applied as short pulses, of 20 seonds 
duration. These pulses of pharmacological agents are indicated by black triangles in Fig. 4a. 
Marked phase increases occur when glutamate (GLUT) is applied alone, at instants indicated 
by black triangles 41, 42, 43 and 44. The concentrations of GLUT, are respectively 15 |iM, 10 

10 jaM, 15 jaM and 20 jiM. At instants indicated by black triangles 45 and 46, 400 |iM of 
kynurenate (KYN), a broad spectrum ionotropic glutamate receptor antagonist, has been 
added to GLUT. No increase of the phase signal occurs in this case, indicating that the 
observed effect of GLUT is inhibited by KYN. The effect of GLUT is also inhibited when 5 
yiM of MK-801, an NMDA specific receptor antagonist, has been added to GLUT in the 

15 physiological solution, at instants indicated by black triangles 47 and 48. As can be seen in 
Fig. 4a, the application of glutamate pulses induced reversible OPL or phase increases lasting 
for several minutes. These OPL increases could be reproduced several times over periods of 
more than three hours. Figure 4a shows also that these glutamate-mediated responses were 
inhibited by kynurenate, a broad-spectrum ionotropic glutamate receptor antagonist, and by 

20 MK-801, a specific NMDA receptor antagonist, thus stressing the receptor specific nature of 
the effect. A concentration-response, curve (Fig. 4c) has been established on the basis of 5 
different experiments. As response amplitudes vary from cell to cell, results were normalized 
with respect to the maximal response which occurred at a concentration of -15 )iM. Curve 
fitting using the Hill equation indicates a highly cooperative behavior (Hill coefficient = 3.75 

25 ± 0.4) with an EC 50 of -7.5 ± 0.23 |iM. Interestingly, these results have been obtained in the 
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presence of glycine (3 ^M) in the perfusion solutions. Without glycine cells started to respond 
at glutamate concentrations higher than 15|aM and for only -25% of cells (vs. ~ 80% in the 
presence of glycine). As glycine is an NMDA receptor co-agonist, this behavior further 
emphasizes the role these receptors in the glutamate-evoked responses. The maximum 

5 amplitude of a typical phase increase is -10 degree representing -5% of the initial cellular 
phase. If we assume once again a constant cellular refractive index (n s = 1.365), this 
corresponds to an increase in the vertical (z) axis of - 0.6 urn. 

Figure 4a shows an example of the evolution of the OPL averaged over a specific area 
of the sample, and gives only a partial idea of the real potential of the method of the present 

10 invention. With DHI, the spatial distribution of OPL can be defined, with a transverse 
resolution of for example 0.3 jam, and OPL time courses, such as presented in Fig. 4a, can be 
obtained for different micro-domains of the cell, or for each pixel of the quantitative phase- 
contrast image. The recorded signal can therefore be followed over time and space, for 
example by constructing animated sequences, for example by using the processing unit 18. 

15 An informative way to characterize the observed effect consists in mapping the spatial 

distributions of OPL variations at different times. For the experiment presented in Fig. 4, this 
representation indicated that the observed cell morphology change is anisotropic. Indeed we 
have observed that OPL measurements presented a marked decrease along the cell 
boundaries, and a global increase over a central region of the soma, on top of heterogeneous 

20 local variations. If we consider an ellipse inside the cell body, the observed shape changes 
consisted in a reduction of its semi-minor axis, which enhances the cell eccentricity. Ten 
minutes after the end of a glutamate pulse application, we observed that small local OPL 
variations persist in some areas, even though the cellular OPL has globally recovered its 
initial level. 
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The high sensitivity of the present invention has revealed that glutamate applied during 
a few seconds at physiological concentrations induces reversible morphological changes 
mediated by NMDA receptors, in microdomains of mouse cortical neurons. These results 
demonstrate that the present invention is a new method that offers unique possibilities for 

5 real-time 3D functional imaging with a sub-cellular level of spatial resolution. 

Morphological neuronal changes (following pharmacological medium applications: 
Glutamate, MK801..) have been monitored during a few hours. To optimize opto-mechanical 
stability during such periods of time the whole holographic microscope may be enclosed in a 
container. Despite this enclosure, different noise sources remain (such as coherent noise), 

10 extremely sensitive to opto-mechanical instabilities, which may perturb the reconstructed 
images. These residual noise sources render the quantitative analysis of neuronal moiphology 
changes very difficult. To reduce the effects of these residual noise sources we have 
developed a device that generates and delivers a controlled energy in the area of the biological 
sample. This controlled energy may be mechanical vibrations and may be transmitted by 

15 perfusion line to the chamber. 

The method of the present invention uses in priority a DHI technique to measure the 
phase of a radiation, and in particular the method described in patent WO 00/20929. 
However, it is evident that all the concepts described in the present invention, and in 
particular the fact that the sample is embedded in a medium, can be extended to other 

20 techniques for phase measurements. In particular, the present invention is also related to 
methods using phase measuring techniques such as in-line digital holography, or Fourier- 
transform digital holography, or digital heterodyne holography, or phase-shifting 
interferometric techniques. Considering that the concept of hologram is a generalization of the 
concept of interferogram, or that the concept of hologram is a particular example of 
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interterogram, any kind ol lnterterometnc tecnniques can oe usea in principle, wunin xne 
context of the present invention, as a mean to measure the phase of an electromagnetic wave. 
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Claims 



1 . An. apparatus for performing digital holographic imaging (DHI) of a sample, said 
5 apparatus comprising : 

- an holographic creation unit, 

- an holographic reconstruction unit, 

- a processing unit, 

- a sample unit, 

10 * characterized by the fact that the sample unit includes a container adapted to 

contain a medium in which a sample is located. 

2. An apparatus according to claim 1 wherein the container is adapted to contain a 
liquid. 

3. An apparatus according to claim 1 or 2 furthermore comprising a medium 

15 controller unit connected to the sample unit, said medium controller unit being 

adapted to modify the medium properties, including the partial or complete 
exchange of medium. 

4. Apparatus according to anyone of the previous claims furthermore comprising 
vibrating means for mechanically vibrating at least a part of the sample unit. 

20 5. Apparatus according to claim 4 wherein said vibrating means are incorporated in 

the medium controller unit. 

6. Medium containing a sample adapted for the apparatus of anyone of the previous 
claims characterized by the fact that it satisfies at least one of the following 
conditions : 

25 - direct effect on the sample, namely by modifying the sample shape and/or its 

dielectric properties, 

- inducing contrast modifications of the holographic images obtained. 

7. Method for using the apparatus of anyone of previous claims 1 to 5 comprising the 
use of a medium, saidmethod being characterized by the fact that the medium is 

30 modified during real-time observations. 
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8. Method according to claim 7 characterized by the fact that the medium is a 
dispersive medium and that the dielectric properties of the medium are changed by 
changing the wavelength of the radiation illuminating the sample. 

9. Method according to anyone of previous claims 7 or 8 wherein the dielectric 
properties of the medium are modified, the corresponding holographic information 
obtained being used to determine the shape and/or the dielectric properties of the 

sample. 

10. Method according to anyone of the previous claims 7 to 9 wherein the orientation 
of the sample is changed during the acquisition of a plurality of holograms and 
wherein the information content of said plurality of holograms is used to build a 
three dimensional computed tomography of the sample. 

11. Method according to anyone of the previous claims 7 to 9 wherein the wavelength 
of the radiation illuminating the sample is changed during the acquisition of a 
plurality of holograms, and wherein the information content of said plurality of 
holograms is used to build a three dimensional computed tomography of the 
sample. 

12. Method according to anyone of the previous claims 7 to 9 wherein both the 
wavelength of the radiation illuminating the sample and the orientation of the 
sample are changed during the acquisition of a plurality of holograms and wherein 
the information content of said plurality of holograms is used to build a three 
dimensional computed tomography of the sample. 
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